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1. Introduction 

Measurements of many environmental variables 
such as wind speeds, rainfall and the concentra- 
tions of atmospheric and aquatic pollutants are 
generally duration-specific; the actual quantities 
measured are averages over a specific time interval 
rather than instantaneous values. Thus concentra- 
tions of ozone are typically measured as parts per 
100 million averaged over an hour, and wind 
speeds are routinely recorded as hourly or daily 
mean speeds. In practice however, the scientists 
wishing to understand the environmental processes 
which lie behind such measurements, and the reg- 
ulatory body which monitors pollution, often wish 
to deal with characteristics measured over some 
other time interval: peak concentrations over a day, 
for example, or high wind gust values, which in 
practice correspond to peak 3-5 s averages. There 
is therefore a need to understand the relationship 
between the statistics of environmental measure- 
ments averaged over different time scales, and to 
relate these statistics to extreme levels of the vari- 
ables. In this paper, results on the asymptotic joint 



distributions of extreme averages over different 
time periods will be treated. These results wttl 
make it possible, for example, to link long historical 
data series containing information about extremes 
of daily rainfall (sometimes extending back to the 
early years of the century) to the shorter series of 
extreme hourly rainfall which have been recorded 
only in the past 20 or 30 years. Thus important his- 
torical information could properly be taken ac- 
count of in the estimation of floods, something 
recognised as higly desirable by hydrologists. An- 
other area of application of the results is in the 
study of the dispersal of airborne pollutants. Here, 
it is known (see, for example Fakrcll and Robins 
[9]) that instruments used to measure the concen- 
tration of pollutants dispersing in a turbulent flow 
cannot resolve the finest scales present in such 
flows. Mea.surements of concentration are there- 
fore invariably obtained only as averages of the 
characteristics of primary physical interest, and so 
a statistical theory which links extremes of averages 
over different ranges would be of great value 
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to scientists working in the area (Mole [10]). One 
related area is the study of joint distributions of av- 
erages and maxima of random sequences. Interest 
in such distributions is motivated by analysis of ex- 
treme winds. For purposes of building design or 
public safety, it is often important to estimate the 
speed of the most extreme wind likely to occur at a 
particular location over a period of years and to do 
this it is natural to apply the methods of extreme 
value theory to data on maximum gusts. The preci- 
sion of estimates obtained may be low due to the 
limited amount of relevant data— often no more 
than 10 or 20 years, so it is desirable to try to im- 
prove the precision by introducing into the estima- 
tion procedure other information relevant to 
extreme winds. It is natural, thefore, to ask whether 
the data on gusts could be augmented by that on 
hourly means, as gusts and means are evidently re- 
lated. One source of guidance here may be provided 
by the limit properties of the joint distribution of 
means and the maxima. 

In Sec. 2, we give a summary of the results that 
exist on the joint limiting forms of sums and maxima 
of stationary sequences and in Sec. 3 we give some 
results on the asymptotic joint distributions of 
extreme averages over different time-periods of 
sequences which have moving average representa- 
tions. Possible solutions for the general stationary 
case are also indicated. 

2. Extremes and Averages 

Let {Xi} be a stationary sequence of random 
variables with marginal distribution function 
P(X^r) = F(jt)andlet 

S„ = S ^'. M„ = max Jr„ , « = I, 2, . . . 

We study here the joint limiting distributions of 

/r i't \ (Sx—tibn M„—d„\ ,,x 

(5„,A/„) = ^^— ,— ^-~J (1) 

as n -» 00 for suitable constants a„ >0,c„ > 0, b„ and 
d„. 

Case 1: Light Tailed Case 

Assume that Vor(Xi) < » and F € £> (A) or F e 
£>{*„), a>2 or F E £►(*»), a>0, where A(.), 
*«(■)' ^"(O ^^^ respectively the Gumbel, Frechet 
and Weibull distributions. 



Then for the associated iid sequence (Chow and 
Teugels [5]) 



or 



(5„,A/,)-''(iV,A), 



where the limit components are independent. 

Can dependence amongst Z, modify this limiting 
independence? As the following theorem (Ander- 
son and Turkman [1]) shows, under quite weak con- 
ditions, dependence does not affect the limiting 
distribution. 

Ttteorem 2.1 

Assume that {Xi} is strong mixing and has posi- 
tive extremal index and for some a„, c„ and d„, 

s„ = 1^ -»" JV(o,i), 



Cn 



where G =A or G =«t'a, for some a>2 or G ='*I^a, 
for a>0. Assume further thatjAi} satisfies the con- 
dition 



lim /:limsup Z>'(an,"-i) = 0, 



where 



ulk 



D'ia„,u„)=2 = E[\exp{ita,r' 2 X,) 



I"! 



-1I[X(A}>«„)1 



(2) 



Then Sn and M„ are asymptotically independent. 
Local dependence condition D'{a„,u„) is quite 
weak and satisfied, for example, by m-dependent 
sequences and by Gaussian sequences with 
summable covariances. 

Case 2: Heavy Tailed Case 
Assume that 

\-F{x)=px-'' L{x), 
F{-x)=qx-'' L{x) 
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where L{x) is a slowly varying function and 
0<a<2. Then for the associated iid sequence 
(Chow and Teugels [5]) 



(S„,M„)^''{U,V), 



where 



Here x(^) denotes the indicator function of the 
event A, Wa(t,p) is the characteristic function of a 
stable law of index a and parameter p, and k is a 
constant depending on a and p. Note that for the 
heavy tailed case, U and V are dependent. Can the 
type of local dependence of the JT-sequence make a 
difference to the Chow-Teugels limit? It can! If 
large values are cancelled by large negative values, 
then sums and maxima can be asymptotically inde- 
pendent. We show this by constructing an example: 
Let {Yi} be a stationary sequence with 

l-Fy(y) = l-y-\ 
e<l,y^l and for v>0 such that t + v< 1, let 



--{-L 



with probability Yi.\ 
with probability 1 - YiZ\ 



Let {£,} be an iid unit mean exponential se- 
quence, T/ = i K, and {5/} iid taking values + 1 and 
- 1 with proijabilities p and q respectively. Then 

Theorem 2.2 

(i) If 0<a!<l, p>0 and conditions D(a„), 
D'{a„) of Davis [6] hold then 






where D =m'm\j: 8/ = I}. 

(ii) If lsa<2, p>0 and conditions D(fl„), 
D'{an) and D'{a„) of Davis [6] hold then 

(S^.H.)-*'' 



/= 1 - 



(iii) Under the conditions of (i) and (ii), 



£:[e-^-x(M„<x)]^'' W^{t,p) <l>„(jt )?-*('■'', 



Then it can be shown that {Xj} is stationary, 1- 
dependentandl-FA'(x) = e(e + v)-'x-%/=>(-x) = 
X"*. Hence the limit distribution of {S^,M„) for the 
associated iid sequence is the Chow-Teugels limit 
with a=e and /» =e(2e+v)"'. The components of 
this limit are dependent. However, (5„, Af„) of the 
dependent (A")} process can be shown to be asymp- 
totically independent due the cancellation of large 
positive values by large negative values values, thus 
showing that local dependence may make a differ- 
ence on the limit distribution. However, if we rule 
out this type of cancellation, then the limit distribu- 
tion is not affected by the dependence in {X}, as 
the following theorem demonstrates (Anderson and 
Turkman [4]). One possible local condition which 
rules out this type of cancellation is Davis' [6] 
T)'(a„) condition, which we assume in the theorem. 
This restrictive technical condition also rules out 
clustering of large and small values above and be- 
low certain thresholds. Types of processes which 
satisfy this condition (and others which we need in 
the theorem) can be seen in Davis [6]. 



where 



/i(r,x)=f [e"<">''--.-i]d(-y-'-). 



and 



5=p 



2-a 



a 



and C is a positive constant. Note that the value of 
the limit does not depend on the dependence struc- 
ture of {Xi}. 

These results seem to be discouraging for statisti- 
cal applications. For example, for sequences with 
finite variance, the independence of (t/, V) does not 
offer a basis for the use of average wind speeds in 
inferences about gusts, contrary to the evidence 
shown in data. This may be due to: 

(i) Time intervals are not long enough in prac- 
tice for asymptotic results to give adequate approx- 
imations, 
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(ii) The correlation structure of the data is not 
well represented by our mixing and local depen- 
dence conditions, 

(iii) Residual seftsonality remains in the data. 
Based on these possible deviations, statistical mod- 
els establishing connection between means and ex- 
treme events are suggested in Anderson and 
Turkman [3]. 



3. Extremes of Averages Over Different 
Time Scales 

The specific problem to be adressed in this sec- 
tion is as follows: Let Jf, represent the instantaneous 
value of the environmental variable at time t, and 
denote by Xj., a moving average of {X,) over the 
range T: 



1 '-' 

Xt,i — '^ 2j X-i' 



(3) 



where 



2 ^(i-.-'X.yi) ' 

k-l 

,.. fl ifj:Gy 
f^^=l0 if not. 



&A 



It is known thatP,, converges weakly as /j-*oo to 
a point process P, a Poisson process with mean 
measure tion R* x{R* -{0})- (Davis and Resnick 
call this Poisson measure, the Poisson random mea- 
sure PRM(^ii.) and consider a more general case 
which involve the left tail of the distribution as 
much as the right tail. Here due to the special sim- 
ple form of the moving averages, we restrict 
ourselves to the space (O.oc).) 

Here 



dfjL^dfXooC'-'' ex(0,«)rfj: 



Hence 



Then we are interested in the paired series {Xt.i, 
Xx.i} for different fixed S and Tand in particular, in 
the joint distributional properties of extremes of the 
^■a\t{Xt.,,Xs.,). 

We will give results only for sequences with the 
heavy tailed distributions. The light tailed case is 
more complicated, since in this case large values of 
the moving averages may occur due to the contribu- 
tion of several relatively large values of the se- 
quence in contrast to the heavy tailed case when 
large values of the moving averages are dominated 
by the largest value of the sequence. The techniques 
to be used to study these questions will be develop- 
ments of those used by Davis and Resnick [8]. 

Suppose that {^/}i*i are iid random variables 
with 



P{X,>x)'-x-'' L{x), 

where L(j:) is slowly vaiying asj:-*», that isA' E D 
(^o). Take constants a„ such that 



nF{X,>a^)'^x~\ 



Consider a point process P„ which puts points at 
(-, (?„■"' X/k), A: = 1, 2, . . . . Hence P„ is the random 
point measure on sets in R* xR 






where {(ujk)} i^ = ] are the points of P. 

Davis and Resnick [8] show that, correspond- 
ingly, for {Xt,.} with same normalization 



2 ^(i.flr'Ar.o (■) -^ Z, Zj CcM.-^i) (■) 

i = l " ' i = l 1=1 



on ;?*x(/?+-{0}), where 



L otherwise 



and (tt,jt) are as above. 

Results for the point process generated simulta- 
neously by {Xt.i} and {Xs,i} processes can be ob- 
tained with a straightforward generalization: 



2 e.i-,..-. 



(.'"n 'Xr,i.a„ 'ATja) 



{•) ^ 2. Z, S'i.;«fi,>«ci)t-) 



(4) 
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on;^*x(/^*-{0})^ where 



c.= |i -1.2.. ...5. 
, otherwise 



The joint limit distribution of Afr.n and Ms.n now 
can easily be calculated since 

P(Mt.^ s a„x, Ms.„ s a„y) 

((0.1]x((-co,x]x(-^.y])')-o) 

^P( E X €fe;V.«^.) ((0,1] X ({- ^,X] X 

(-co,y])^=0) 
=P(2 €(...r-.;,.5-„,((0,llx((-«>,J:]x 

(--,>'))') =o) 

= exp[-n((0.11x(r;cAii>',"))] 
= exp[- ?(jr)dx] 

JTlASi' 



'TtASy 



'-{ 



exp [ -p (7x A 5k)-°] if Tj: a S>' > 
if not 



(5) 



Special Cases 
[1] 5 = 7". Then 

/'(a,r'Mr.. :Sx, a,,' Ms,„ <y)=P(a,r'A/7:« ^xr\y ), 

so from theorem 3.1 of Davis and Resnick [8], we 
should have 



Hence the result Eq. (5) is consistent with the 
existing results on this special case. 

f2] S = \<T, ^\Ms.n =max,sr. A^-, and the result 
Eq. (5) says that 

P(a„-'MT,n ^x,a„-'Ms,„ <y)^e-''^'^^>'^-\ 

Note that the above distribution is the joint limit 
distribution of the maximum of the {Xi} process and 
the maximum of a moving average of it. 

The above set up is a very simple case. The imme- 
diate question is: what if the {Xt} are themselves a 
dependent sequence? 

We can get a partial answer by taking {X,} to be 
itself a moving average: 

EC 

7 = 1 

say, where Z, satisfy the same conditions as the Xi 
and fl;SO (for simplicity) 
Thus 



and 



Yr..={X}*[^, 
X = {Z}*{a}. 



Hence 

say, where 

Note that for any 1-1,2, ..., 

^ TAI-l 



As before 



Pion ' Mr.n ^x, a„"' \fs.n ^y) 
-*HS 2 C(.w..*..;,<i,)((0,l]x 

((-=c,j:]x(-^,y])') = 0), 



(6) 
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where 



Note that 






&iuj.^.ni>mnxii~- « .^1 X ( - «.;'])') =o}s 



for every i for which one of the compKJnents is non- 
zero. This is the case iff jkd^^x andj**/* <y, where 
d* - max d; and d* = max di. Hence the limit of Eq. 
(6) is given by 



p(2^AjtAJt.^)=o) 






10 otherwise 



(V) 



Hence we see that this kind of dependence in the 
underlying process does not change the form of the 
limit distribution. It would be interesting to obtain 
similar results for general stationary sequences. 
This could be done by using characterization of the 
limit point processes for the sequence of point pro- 
cesses with points {t/n,a^' X,,t = 1,2,... , n) given 
by Davis and Hsing [7] when {X,} is a stationary se- 
quence with regularly varying tails, in their paper, 
Davis and Hsing show that when {X,} satisfies a 
proper mixing condition then 






converges to a point process N of the form 

to as 

/-I v-l 



Zr Zr e(5,.p,) 

/-I ;-l 

is PRM(v) with 

dv = dtxyax~''~^dx,x>0, 

(This is slightly weaker form of Davis-Hsing limit, 
since we consider the convergence only on R"^ -{0} 
not involving the left tail.) They also show that un- 
der the proper mixing condition, the convergence of 
N„ to N is equivalent to: 

For r«-*«', r„/«-*0 as n-*a>, and k, =^[n/r„]. 



lim k„ P( max Xi>a„x) = -yx '", 



and 



^(2 ^Xiimatti^r^xj G./max Xj>a„x) 



ISjSfn 



/=1 



>Q() 



Here Q(.) is the distribution of the iid point pro- 
cesses 






From this basic result it may be possible to obtain 
the limiting form of 



^"-t^(i;,a„"Xr.,.a„-Xs.,)(-)' 



where 



and 






1 ^ 

Xs,i ="p 2j -*'-;■- 



where 
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One would expect that if the above sequence of 
point processes converges, then the limit point pro- 
cess should be of the form 



n; 



Zj 2j Z, ^lS„ciP,Q,^elF,aj), 



where 



t n 



otherwise 
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1 other 



= 1,2. ..,S 
otherwise 



Straightforward adpstment of Resnick's [11] 
arguments which involve consecutive application of 
the continuous mapping theorem is not possible, 
since these arguments use the convergence of 

nP{fl„-' Xi>S, a„-' X2>S)-^(i,S>0 

which is clearly not satisfied by most sequences with 
strong local dependence. 
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